A CLASS OF MATRIX- VALUED POLYNOMIALS GENERALIZING 

JACOBI POLYNOMIALS 
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^^ , Abstract. A hierarchy of matrix-valued polynomials which generalize the Jacobi poly- 

fT^ ' nomials is found. Defined by a Rodrigues formula, they are also products of a sequence 

OA , of differential operators. Each class of polynomials is complete, satisfies a two-step recur- 

rence relation, integral inter-relations, and quasi-orthogonality relations. 
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r^ ! 1- Motivation 

-)— > . 

J2 ' The understanding of matrix-valued orthogonal polynomials has advanced greatly in 

recent years, and a wealth of references can be found in the recent work of Barry Simon 

The polynomials introduced in the present paper generalize the Jacobi polynomials. A 
K* I different generalization, yielding an orthogonal class, is due to Griinbaum [2]. 

The polynomials defined here were not sought-for. Rather they appear naturally in the 
f — ! . forefront of a different problem, the study of the possibility to linearize (by an analytic 

cn ' change of variables) a differential equation whose linear part has only regular singular 

'!;:j- ■ points. In a neighborhood of one such point equations are (generically) linearizable [3]. 

^^ I But this is no longer the case if the domain studied contains two such singularities. For 

2^ ■ example: 









(1) -^ = Mu -h f(x,u) withM = A+ B 

ax X — 1 a; -I- 1 



H ' (f collects all the nonlinear terms in u) is not necessarily equivalent to its linear part 

^ '. 

, , dw 

(2 -— = Mw 

ax 

for X in a domain in the complex plane including both singular points ±1. 

It turns out, however that for any nonlinear term f (x, u) there exists a unique 4>{u) so 
that the equation with the "corrected" nonlinear part f(x,u) — (/'(u) is linearizable [4j. 

Besides its clear intrinsic interest, the problem of detecting linearizable equations is 
important also since linearizability and integrability turn out to be intimately connected 
0. 



series m 
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Looking for an analytic change of variables y = w + ii.{x, w) that transforms ([2]) into 
([1]) one is lead to the study of the homological equation which boils down to equations of 
the form 

(3) (9^h + (iwhMw-Mh = g(x,w)/(x2 -1) 

where now all functions can be assumed homogeneous polynomials in w, of degree, say, n. 

The main question is: for which functions g does equation ^ possess a solution h which 
is analytic in a domain containing both singularities x = ±1? 

In the scalar case, when w = 'w£C,h = h,g = g(^C, and A = a, B = b & C, it turns 
out that: equation ([3]) has an analytic solution if and only if the first coefficient in the Jacobi 
I pU" ;« A" j )[ of 5 is zero. Moreover, Jacobi series expansions are 

canonical to approaching the problem, since for g{x) = -P^V\ (x) the unique 

analytic solution of dSD is /i = ^p««-i)«'("-i)^). 

In the linearization problem mentioned the double sequence of Jacobi polynomials p^ "■"'"■ -' 
for A; > 0, n > 2 is interconnected. Numerical expansions in Jacobi series yield (what appear 
to be) rapidly convergent seriesU 

The natural question was whether the polynomial structure found in the one-dimensional 
case survives in more dimensions. 

The answer turns out affirmative - see (j28p . (j29|) . These multidimensional polynomials 
share many properties common to the usual classes of orthogonal polynomials: each class 
is complete, satisfies a two-step recurrence relation, quasi-orthogonality relations exist, and 
there are integral inter-relations. It is the author's belief that more that what is presented 
here does hold. 

2. Main Results 

2.1. Notations. Consider a Fuchsian differential equation in C^ with three singularities 
in the extended complex domain C: 

(4) y' = My where M = -^ A + ^— B , (xeC, ye&) 

X — 1 X + 1 

where A and B are d x d matrices so that 

(i) A + B is diagonalizable, and 

(ii) A + B satisfies the following nonresonance condition: the eigenvalues Ai, . . . , A^ of 

the matrix A + B satisfy: 

(5) n • A - Aj Z_ for ah n € N'^ , j = l,...,d 

These assumptions are, of course, satisfied by generic matrices. 
Denote for simplicity: 

(6) Mi = A + B, M2 = A-B, Q = x'^-1 



However, the proof of convergence of the series obtained, which did not use the Jacobi structure, turned 
out to be extremely delicate: a steepest descent method was used and small denominators had to be dealt 
with [6]. 
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SO that QM = xMi + M2. 

It will be assumed, without loss of generality, that the matrix Mi is diagonal. 
Let y be a fundamental matrix of solutions of @ , therefore 

Y' = MY and {Y'^)' = -y^M 

Denote by Vn the space of vector- valued polynomials in vi^ = (wi, . . . , Wd), homogeneous, 
degree n: 



Vn = { q ; q(w) = Yl qmw™, q„, G C" 

mGN'^,|m|=?i 

where |m| = mi + . . . + rrid, w™ = w'[^^ . . . w^"^ , N = {0, 1, 2, . . .}. 



2.2. The class of polynomials: definition. The polynomial in x degree k, Pfc(x), is, 
for every x, the linear operator on Vn defined by 



(7) 



Pfc(x)q(w) = y(x) 






t=x 



Q(t)^y-i(t)q(y(t)y-i(x)w) 



Formula ([7]) can be compactly written as follows. Denoting by W{x) the linear operator 
on Vn- 

W(x)q,{w) = Y{x)-^ci{Y{x)w) 

(which gives q in new, x-dependent coordinates in C^) the definition ([T]) is 

d^ - 



(8) 



Pfc(x) = W{x) 



-1 



dx^ 



Qix^Wix) 



which is a Rodrigues formula. 

The fact that Pfc(a^) are polynomials in x, degree k, is shown in Proposition[TJ 
Denote by Di , D2 the following linear operators on Vn '■ 

(9) Aq(w) = dq(w) M^w - Af,q(w), i = 1, 2 

where Mi are defined in ([6]). 

Proposition 1. The operators Pfe(x) are polynomials in x degree k. 
Moreover, they can be written as the composition 



(10) 
wher^ 

(11) 



A1A2 ...Ak 



Aj = x (2j + Di) + D2 + Qd^ 



Note that the action of first operator Ak when apphed to q(w) (which does not depend on a;) equals 
[x{2j + D^) + D2Mw). 
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Proof. 

Note that calculating the first derivative in ([7]) we obtain 

d 



dt 



QityY{t)-'(i{Yit)u) = Q{ty-'Y{t)-' {2kt + tDi + D2) q(y(t)u 



Q{t)''-^Y{ty^r{t, Y{t)u) where r(t, w) = Aq(w) 



Then 



dt'^ 



Q{tfY{£)-\{Y{t)M) 



d_ 

IE 



Q{tf-'Y{t)-'v{t,Y{t)M 



= Q{tf-'^Y{t)-^ [2{k - l)t + tDi + D2 + Qdt] r(t, Y{t)u) 
= Q{t)'''^Yit)-Uk-ir{t,Y{t)n) = Q{tf-^Y{t)-^ Au-iAMY{t)vi) 

The k derivatives in ([7]) can be calculated in this way recursively, yielding (jlOp and 
completing the proof of Proposition [TJ D 

Remark 1. The same definition ([7]) can be used to define Y'j^{x) as linear operators on 
the space of C^-valued formal series in w, or on the space of convergent such series. 

2.3. The polynomials ([7]) generalize the Jacobi polynomials. Consider the one- 
dimensional case: d = \. Equation (j3|) and a fundamental solution are, in this scalar 
case, 

a b 



y 



+ 



y 



and 



y(x) 



l)'^(x + l) 



X — 1 X + 1 ^ 

The homogeneous polynomials degree n are just multiples of w": q{w) = cw"' and formula 
dZ]) gives 

1 d'' - 



Pkix) 



y{x)"' ^ dx^ 



Q{xry{x[ 



,n-l 



which is the Rodrigues formula for the Jacobi polynomial -P^ {x) (up to a 

multiplicative factor) . 

Remark 2. In the one-dimensional case the operators D12 of ([H) are multiplication 
by (n — l)(a -|- b), and (n — l)(a — b) respectively, and formula (fTO]l gives the following 
representation for Jacobi polynomials (up to a numerical factor): 



Pt'\x) 



A1A2 ...Akl 



(12) 

where Aj = (2j + a + P)x + {a - P) + Qd^, 

Formula (1121) provides the following elegant way to deduce the second order operator for 
which Jacobi polynomials are eigenvalues: noting that 

d^Aj = {a + P + 2j) + Aj+A 
we obtain 

AAP^^^ = Aid^AiA2 ...Ak = {a + P + 2)P,("'^) + ^i^aS.^s • • • A 
= . . . = [(a + /3 + 2) + . . . + (a + /? + 2k)] PJf 



k{a + p + k + l)Pt''''^ 
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2.4. The multidimensional commutative case. If the matrices A and B are simulta- 
neously diagonalizable: 



A = diag [aj] •=! ... d ™'i ^ = diag [6i]j=i 



then the matrix Y is 



Y = dias 



(x-lf'(x + l)^ 



i = 1, . . . , d 



and the action of the operators ([7]) on the canonical basis of Vn- ^m,j = 'w^'ej , j = 1, . . . ,d, 

Iml = n is 



(13) 

where 
(14) 



Pk(x)h- 



1 



d" 



mj 



VFinj(a^) dx'^ 



Q{x)''W^,,{x) 



-"mj 



W^m,,(x) 



(x-1) 



(x + 1 



imb— 6, 



therefore Pa;(x) are diagonal operators, with entries Jacobi polynomials degree k. 

This is a family of matrix- valued orthogonal polynomials with respect to the inner prod- 
uct 

<P;,Pfc>:= / Ti {Pi{x)Wix)Pkix)} dx 

\m\=n,j=l,...,d 

in the cases when the integrals exist o 

2.5. The non-commutative linear case. For n = 1, the homogeneous polynomials 
degree one are linear transformations: q('w) = qw where g is a matrix. Formula ([T]) is, in 
this case, 



(15) 

which implies 

(16) 



d 



Pk{x)q = Yi^)^ [Q{x)''Y-\x)qY{x)\ Y-'{x) 



Tr [Fkix)q] = Pi°'°^(x)Trg 



where P^ ' (x) = Pk{x) are the Legendre polynomials. 
2.6. Properties of the generalized polynomials ([?!) . 

2.6.1. Rodrigues formula. The definition ([8|) of Pfc(2;) is a Rodrigues formula. 



Appropriate conditions relating a, b, n, d will ensure existence of the integrals. 
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2.6.2. A two-step recurrence relation. 

Proposition 2. The operators Pfc(x) defined by ^ satisfy the following two-step recur- 
rence relation 

(17) xPk{x) = Pfe+i(x)afc + Pfc(2;)/3fc + Pfe-i(x)7fc 
where ak,/3k,^k o.re linear operators on Vn, and in particular 

(18) Ofc = {Di + 2k + ly^Di + 2k + 2)-\Di + k + 1) 

Proof of Proposition^ 
Note the identities 

(19) Aj{xy) = xAj{r) + QT 
and 

(20) Q Aj{r) = Aj-i{Q r) for any r = r(a;, •) G Vn 
which by iteration give 

A1A2 . ■ ■ Ak{x q) = A1A2 . ■ ■ Ak-i {x Aq) + A1A2 ■ ■ . Ak-i {Q q) 

(21) = ... =x A1A2 ...AkCi + k A1A2 . . . Ak~i {Q q) 
Using PU]) and ([21]) relation p!7|) follows if we have 

AkX -kQ = AkAk+iak + A/3fc + 7fe 

which expanded yields an identity of quadratic polynomials in x, and by identifying the 
coefficients we obtain the following equations for Ok, I3k,7k'- 

(22) k + l + Di = {2k + l + Di){2k + 2 + Di)ak 

(23) D2 = [{4:k + 2)D2 + D1D2 + D2D1] ak + {2k + Di) f3k 

(24) k-l=[Dl-{2k + 2 + Di)] ak + 1^2^ + Ik 



The system (l22]) - (i23]) has a unique solution ak,Pk,1k due to the nonresonance condition 

dSJ. □ 

2.6.3. The dominant coefficient. 

Proposition 3. The coefficient of x in Fk{x) is 

(25) {Di + k + l){Di + k + 2) ...{Di + 2k) 

which is invertible under the assumptions of ^2.1[ 
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Proof. Retaining only the dominant coefficients in the representation (jlOp . (jlip we have: 

Pk{x) = [x(2 + L»i) + x29^] ... [x{2k-2 + Di) + x^d^] [xi2k + Di)] + O(x^-i) 
which yields 

Pkix) = iDi + k + l){Di + k + 2) ... {Di + 2k) x^ + O(x^-i) 



giving (j25]) . 

Under the assumptions and notations of (i), (ii) of §2.11 we can assume Mi diagonal, 
and the operator Di has the eigenvector /values: 

Dihmj = (m • A - Aj) bmj 

Then Di + j is invertible for all j € Z+ due to the assumption (ii) , and therefore ([25]) is 
invertible. D 

2.6.4. Completeness. The set Pfc(x), A; = 0, 1, 2, . . . is complete in Vn[x] in the following 
sense: 

Proposition 4. For any f = f(a;, w) E Pnl^;] polynomial, homogeneous degree n in w and 
degree k in x there exist qo, . . . q^ € Pn so that 

k 
(26) f(x,w) = ^P,(x)q,(w) 

j=o 



and the representation I126\) is unique. 



Proof. The decomposition ([26]) follows easily by induction on k, relying on the fact that 
the dominant coefficients (I25D are invertible. D 



2.6.5. Integral inter-relation. Consider the following variation of the Fuchsian equation ^: 

2x 

and let y be a corresponding fundamental matrix: Y' = MY. Let Pk denote the polyno- 
mials associated by ([?]). 

We have the following analogue of an integral relation for Jacobi polynomials (see [^ 
§22.13.1) 

(28) PA:(x)q(w) = Q{xy^Y{x) J'' Y {ty^P k+i{t)q_ {Y{t)Y{x)-^w) dt 

or, in differential form, 

(29) 5^Pfcq + d^(Pfcq) Mw - MP^q = Q'^Pfe+iq 
or, in operator notation, 

(30) {xDi + D2 + Qd^) Pfe = Pfe+i 
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This can be easily seen due to the fact that P^ also satisfy 

Pfc+i(x)q(w) =Qy(x)^^^ ^^^ [Qit)''Y-\t)ci{Yit)Y-Hx)w) 

(which follows by a short calculation using the fact that q(w) is a homogeneous polynomial 
in w, degree n.) 

2.6.6. Orthogonality. The following quasi-orthogonality relations hold: 

(31) / Pj{x)W{x)Fk{x)dx = ioT j <k 
and 

(32) / W{x)Pj{x)Pk{x)dx = for j > k 

provided that the integrals exist. 
Proof. For any index k we have 

d^ r 

[w(x)Pfc(x)]q(w) = y(x)-ipfc(x)q(y(x)w) = ^ [q(x)'= y-i(x) q (y(x)w) 

and (|3ip . (|32p follow using integration by parts (which holds for matrix multiplication). D 
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